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Abstract 

This paper presents a quadratic approximation for the optimal power flow in power distri¬ 
butions systems. The proposed approach is based on a linearized load flow which is valid for 
power distribution systems including three-phase unbalanced operation. The main feature of 
the methodology is its simplicity. The accuracy of the proposed approximation is compared to 
the non-linear/non-convex formulation of the optimal power flow using different optimization 
solvers. The studies indicate the proposed approximation provides a very accurate solution for 
systems with a good voltage profile. Results over a set of 1000 randomly generated test power 
distribution systems demonstrate this solution can be considered for practical purposes in most of 
the cases. An analytical solution for the unconstrained problem is also developed. This solution 
can be used as an initialization point for a more precise formulation of the problem. 

Keywords: Optimal power flow, convex optimization, quadratic optimization, approximation 
models for power system analysis 


1. Introduction 


Optimal power flow (OFF) is a classic problem for transmission system operation which has 
been extensively studied in the scientific literature Hi] . The increasing penetration of renewable 
energies and the possibilities offered Iw communications in future smart-mds allow the use of 
OFF in power distribution systems llslQl and especially in micro-grids 13,Si- 

OFF is a challenging problem due to the high number of nonconvex constraints. Newton- 
Raphson, descendent gradient and interior points methods are traditionally employed to obtain 
suboptimal solutions uJ]. These methods allow decoupled formulations in the context of trans¬ 
mission networks, since nodal voltages are usually close to IZO and reactance/resistance ratio of 
transmission lines is frequently high. A good quality initial solution as well as a simple modeling 
which allows fast calculation of derivatives, are key features for a fast and accurate solution of 
the problem ||3l. 

In power distribution systems, it is a further challenging problem due to the unbalanced 
operation and low reactance/resistance ratio of distribution lines. Therefore, heuristic algorithms 
based on artificial intelligence have been proposed to find good solutions to the problem Hil^. 
Evolutionary algorithms ifToi [Til and particle swarm optimization 112 ] have demonstrated to 
be efficient approaches for the problem. This type of algorithms allow an accurate modeling 
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of the system by including constrains otherwise very difficult to consider. However, heuristic 
algorithms do not guarantee optimality and can be computationally cumbersome for real time 
operation. 

Another approach for the problem in both, transmission and distribution networks, is the 
use of relaxations and simplifications in order to ’’convexify” the problem IB [USE Q- 
Semidefinite programming is one of the most promising modeling techniques for this propose 
118. B- The main advantage to reformulate a problem as a convex optimization problem is 
the capability to find global optimal solutions in an efficient way IBl . In addition, a convex 
formulation allows in some cases, the use of distributed methods. This is a key feature for future 
smart-grids. 

The difficulty of the OFF lies in the non-convex nature of the load flow equations rather than 
in the number of variables. Different convex approximation have been proposed in the literature 
to address this problem. For example, in [210 ^ curve-fitting technique was used in order to 
linearize voltage-dependent load models. Other analytical approaches were presented in [221] 
and IBl . 

This paper introduces a quadratic convex approximation for the OFF in power distribution 
systems. This approximation is based on the linear formulation of the power flow presented in 
12311 . Different consideration are made ending at a non-iterative analytical solution for the relaxed 
problem. Both, the quadratic convex model and the analytical relaxed model are extensible to 
three-phase unbalanced distribution systems. These results have many potential applications 
including: 


• As initial point for other non-linear or heuristic algorithms 

• Asa practical solution in systems in which a close-to-the-optimal solution is acceptable 

• In markets regulation where a convex formulation is desired (i.e the solution is unique and 
do not depend on the used algorithm). 

• In real time operation were a fast solution is required 

• As part of other algorithm that requires to call many times an OFF as a sub-routine. 

• As sensitivity analysis for power distribution systems 


Unlike conventional formulations, the proposed approximation uses complex voltages as 
state variables re pres ented in rectangular form. Although similar formulations have been pro¬ 
posed before |24, 25, B B] the results presented here are different in three main aspects: First, 
the proposed formulation seeks an approximated model rather than an efficient implementation 
of a conventional Newton-based algorithm. This approximation has theoretical and practical ap¬ 
plications from the power engineering stand point. Second, modeling and linearization is made 
entirely on complex variable before split in real and imaginary part for the optimization process. 
Off course, it might be possible first split and then linearize, but modeling in complex variables 
allows a straightforward extension to the three-phase case and inclusion of complex constrains. 
Third, an non-iterative solution is found for the relaxed case. Due to the non-convex characteris¬ 
tic of the problem, even the unconstrained case is difficult and can lead to local optimums IB. 
Therefore, a global non-iterative solution is useful even for initialization purposes 12811 . 

The remainder of this paper is organized as follows. Section II reviews the non-linear non- 
convex formulation of the OFF and analyzes the advantages of a rectangular formulation in 
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power distribution systems. Section III presents the quadratic convex approximation of the OPF. 
In Section IV, further approximation are considered to end at a non-iterative analytical solution 
for the relaxed problem. In Section IV, the methodology is extended to three-phase unbalanced 
systems. Finally, Section V presents simulation results performed over an extensive set of test 
systems before Section VI concludes. 

2. Formulation of the OPF for power distribution systems 

Different formulations for the OPF have been proposed in the scientific literature as a result 
of contributions from many researchers in this area. Two main formulations can be considered 
namely Polar-OPF and Rectangular-OPF, according to the representation of the state variables. 
Both formulations are equivalent. In the first case, decision variables are active and reactive 
power of distributed generators while voltages are state variables represented in polar form. In the 
second case, decision variables are currents injected by generators and state variables are voltages 
(both represented in rectangular form). Different objective functions can be considered including 
minimal generation costs, maximum market surplus and minimum losses, among others. In this 
paper, the minimum losses OPF is considered although the methodology can be extended for 
other objective functions. 

Rectangular-OPF is less common in the literature than Polar-OPF |[ll]. However, it has some 
advantages in power distribution systems, especially in those cases where distributed generators 
are operated at constant power factor. In this formulation, voltages and currents are represented 
in rectangular form ((vr, v,)(!V, /,)) as given in Eqs. ([T]i to (|9]l 
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where subscripts r, i represent real and imaginary part, and subscripts (k), (m) represent nodes 
(with (0) the slack node). Moreover, g(icm), b{icm) represent real and imaginary components of the 
nodal admittance matrix respectively, p and q are the nodal active and reactive power, Pg, Qg are 
the active and reactive power delivered by distributed generators, and finally a is an exponent 
that represent the ZIP model of each load (i.g 0 for constant power loads, 1 for constant current 
and 2 for constant impedance). 

The main advantage of the Rectangular-OPF is that coupling between nodes are represented 
by a linear equation while the non-linear/non-convex equations are isolated to each bus. The main 
source of non-linear equations are constant power loads. Distributed generators can be consid¬ 
ered as PQ buses for mathematical optimization modeling even if operated at constant voltage. 
Set point of the voltage can be calculated after the optimization is performed. On the other hand, 
most of the loads in distribution systems require a model which considers constant power, con¬ 
stant current and constant impedance loads (i.e the ZIP model). Therefore, a linearization of the 
constant power loads is required in order to obtain a convex approximation. 


3. Quadratic approximation 


This section presents the quadratic approximation from the Rectangular-OPF. The key step 
in this development is the linearization of the load flow equations which was first presented in 
I 23n . For the sake of completeness it is briefly presented below. 

Let us consider a power distribution system whose topology is described by the nodal admit¬ 
tance matrix as follows: 


/ lo \ / Yoo Yon W Vq \ 

\ j \ Yno Ynn I \Yn I 


( 10 ) 


where 0 is the substation node (slack) and N — {1,2,...«} are the remained nodes. Along 
this section, a blackboard bold variable represents a complex matrix or vector while an unbold 
variable with sub index r or i represents its real or imaginary part respectively. Slack node is 
assumed at 0° hence Vo = Vo 

Current and voltage in each node k e V are related to the ZIP model of the loads by Eq. (fTTb 


= ^ + + (11) 

Notice that component P of the ZIP model is the only non-linear equation. An approximation 
of this term was proposed in lEill . based on the Laurent series expansion lE^ of the function 
F(V) = 1/V within a closed region ||1 - V|| < a < 1 as follows: 


^-S;,-(2-V:) (12) 

This approximation can be also obtained by spliting Eq. (fTTT i in real and imaginary part, 
and expanding in a Taylor series around V = 1 . However, it is straightforward to maintain the 
formulation on the complex set. It allows a linear and compact equation for the load flow which 
is analog to the DC power flow for power transmission networks as follows: 
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with 


a + i-v;, + c-Va, = o 


(13) 


A = Y^oVo-2-S^-S;^ (14) 

B = diag(S*p^) (15) 

-diagiS*z^) (16) 

It is important to notice that not supposition is made related to whether or not the system 
is radial or balanced. Unlike the DC formulation of the power flow in transmission networks, 
nodal voltages are not fixed to 1 pu but obtained by the power flow itself. Thus, PV nodes are 
not considered. Nevertheless, PV nodes, if any, are less common in power distribution systems 
than in transmission networks. In addition, distributed generation can be considered as PQ nodes 
during the optimization process regardless of the type of control in real time operation. 

The exactitude of the linear formulation of the power flow and its extension to three-phase 
systems was discussed in lE^ as function of the maximum voltage drop (6max)- 

l|l-V||<d,„„, (17) 

The proposed approximation is valid for power distribution systems in which Eq. (ini) is 
fulhlled. For example, a dmax - 0.3 (i.e 0.7 < ||V|| < 1.3) results in a maximum expected error of 
10% which decreases as ||V|| —> 1. In addition, angle for the slack node must be 0°. Otherwise, 
a rotation is required as presented in Section|4]for three-phase systems. 

Let us include the distributed generators in Eq. O as follows: 


A + M-Y*f^ + C-YN + D-Sl;=Q (18) 

where D - {dic,m)i^xc ^ matrix such that dkm = 1 if the generator m is connected to the 
node k, and 0 otherwise. This affine equation replaces the non-linear expressions given by Eqs. 
|2]to|5] The additional convex constrains required to complete the optimization model are given 
inEq. ( fT9l i: 


Minimize - Real -i- 2 • V^Yato • Vo) 

subject to A H- B • VT + C ■ + D - -0 (19) 

l|l-V||<d„„, 

IIScll ^ Scimax) 

After separating in real and imaginary part, all equality constrains are affine functions whereas 
all inequalities are convex function. Symbol ||.|| could represent absolute-value norm or Eu¬ 
clidean norm of the complex number. Both cases result in a convex optimization problem. How¬ 
ever, the hrst case is easier to solve since it is a quadratic/convex optimization problem with 
linear constrains. 

The optimization model given by Eq. (fT9l l is compact and clear due to the representation 
in complex variables. This is important for implementation purposes in a software that allows 
complex variables (i.g Matlab/Octave). Objective function cannot be complex on account of 
it is not an ordered field. However, feasible region and decision variables can be efficiently 
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represented on the complex field. Nevertheless, a more conventional representation of this model 
is obtained by separating in real and imaginary part starting from Eq. (HI as follows: 


I D Sr\ /Ar\ /Br + Cr Bi — Ci 
\-D Sij ^ A, ; \ Bi + Ci -Br + Cr 



( 20 ) 


Notice this expression allows to obtain Vr and Vi as function of Sr, Si (ie. the active and 
reactive power delivered by distributed generators) since remaining terms are all constant. Let 
us define a constant matrix Mn^n as follows: 


M = 


Br + Cr Bi — Ci \ 

Bi + Ci ~Br + Cr j 


then (I2OI 1 can be re-written as 


( 21 ) 



D Sr ] 
-D Si I 


( 22 ) 


let us define U - Ur + j-Ui as the vector voltages of the system without distributed generation 
(i.e Sg = 0). 



and W 2 nx 2 n a real matrix given by 


(23) 


/ Wri 

\ Wir Wii 




(24) 


These definitions allow to represent nodal voltages as function of the generated power. 


I Vr \^l Ur \ ( Wrr Wri W \ 

E,- j t/, j 1T,> j 


(25) 


Nodal voltages V are equal to U when generated power Sc is zero, which correspond with the 
definition of U. At this point Eq. (l25T l is just another representation of Eq. (fTSl lthat is replaced 
on the optimization model ( fT9l l as follows: 


Minimize Vj ■ Gn ■ Vr+ Vj ■ Gn ■ Vi+ 2 ■ Vj ■ Go ■ Vo 
subject to 

Vr^ Ur + Wrr-Sr + Wri-Si 

Vi^Ui + Wir-Sr + Wii-Si (26) 

(1 - Vr(,)f + < dLx 

S r(/c) ^ r(max) 

S i(lc) ^ S i(max) 

where 


Gn — Real(Y nn) (27) 

Go = Real{'¥]^ q) (28) 
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The model is convex and hence, it can be solved easily no mater the size of the system. 
Gn is a square matrix but Ms not. Therefore, Vr - V, = 0 are not feasible solutions. In 
addition, if maximum drop constraint (Eq. [TTI i is relaxed, then the model becomes in a simple 
quadratic optimization with linear constraints which can be solved by well known toolboxes 
such as quadprog in matlab. Nevertheless, Eq[T7]must be checked after the optimization since 
the exactitude of the load flow depends highly of the this constraint. 

3.1. Analytical relaxed model 

Inequality constraints in Eq. |26]can be relaxed in order to obtain an analytical relaxed model. 
This model can be used for initialization of the non-linear optimal power flow as will be presented 
in Section|5] 

Let us include load flow equations into the objective function as follows; 


Minimize (Ur + Mr • Sr + Wri ■ Sif ■ Gn ■ (Ur + Mr • Sr + Wri ■ Si) + 

(Ui + Wir ■ Sr + Wii ■ Si)^ ■ Gn ■ (Ui + Wir ■ Sr + Mv 'Si) + (29) 

2(Ur + Wrr-Sr + Wri ' 5,)^ ' Gq ■ Vq 


which can be re-written as 


with 


Minimize -S^ ■ H ■ S + ■ S 


H = 


Hr, 


Hi, 


Hri 


Hi 


(30) 


(31) 


Hrr = M^r ' Gn ' M^^r + ■ Gn ' (32) 

Hri = M^r ■ Gn ■ M^ + K-Gn-WI (33) 

Hir = wl ■ Gn ■ M^r + wI-Gn-WI (34) 

Hu - Wji • ■ Mv + Gw • Mt (35) 


and F - (Fr, FiY is given by 


Fr = M^r ■ Gn ■ uj + wJ'^-Gn- uJ + -Gq-Vq (36) 

Fi = ■ Gn ■ Uj + M^ ■ Uj + ■ Go ■ Vo (37) 

where angle of the slack node is assumed to be equal to zero (i.e y,o = 0). The optimal 
solution of the relaxed optimization problem (i.e analytical relaxed model) is given by Eq. (l38l l. 

S = -i/-' ■ F (38) 

This analytical solution has many applications. In a general case, it can be used as an initial 
solution for a non-linear optimization methodology or a sensitivity equation for other optimiza¬ 
tion problems. 
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4. Extension to three-phase systems 

The proposed approximations can be extended to three-phase unbalanced distribution sys¬ 
tems with some additional considerations. Nodal admittance matrix must consider a three-phase 
modeling with a three-phase Slack bus. As a consequence, Vo is not a constant but a 3 x 1 vector. 
In addition, Y-connected as well as A-connected loads must be considered. Per unit represen¬ 
tation is not an advantage in this circumstances. Thus, a new constant rf - llv„om is defined in 
order to permit the linearization around 1. 

The basic structure of Eq. [T8]remains but now the matrices are redefined as follows: 


A = Yno -Vq-It] ■ -jj ■ o Tat (39) 

B = 7?^ • ■ diagiSpf^ o T^) ■ (40) 

(41) 

D = 77 ■ ■ diag{TN) ■ D (42) 


where T(Arxi) is a complex vector whose elements correspond to each node such that Tk - 
with (f) e {0,-27r/3,27r/3} according to the sequence. Symbol (o) is the Hadamard product 
(i.g the matlab element-wise multiplication) in order differentiate from the conventional matrix 
multiplication. Moreover, 7 is a real matrix which indicates the type of connection of each load: 
for Y-connected loads, it is an identity matrix while for A-connected loads it is a matrix that 
converts phase voltage into line voltages. 

On the other hand, D has the same definition as in the single-phase case but a new complex 
matrix D requires to be defined as Eq. (l42l) . Additional constrains related to the three-phase 
modeling could be included directly on D. for example, some generators could be inappropri¬ 
ate for unbalanced generation (i.g renewable resources integrated by three-legs voltage-source- 
converters). In those cases, generated power on each phase must be equal despite the unbalanced 
operation of the rest of the system. This is represented by Eq. dmi. 

(43) 

where k, m, n are different phases of the same three-phase bus. This constraint can be included 
directly in D. To take a simple example, consider a three-phase distribution system with 3 nodes 
and 2 distributed generators connected to node 2 and 3 as depicted in EiglT] 



Figure 1: Three-phase distribution system with disti'ibuted generators 
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Generator connected to node 2 can generate three-phase balanced power while generator 
connected to node 3 can generate power independently to each phase. The matrix D for this 
system is given by Eq. (l44l) where nodes are ordered by phase (i.e fist the phase a for nodes 
1,2,3, then phase b and finally phase c). 


0 

0 

0 

0 

1/3 

0 

0 

0 

0 
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0 

0 
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0 

1 

0 

0 
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0 

0 
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0 

0 

0 

0 

0 
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If generators is and loads are A-connected, then Jnn is given by 


1 

0 

0 

-1 

0 

0 

0 

0 

0 f 

0 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

-1 

-1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

1 j 


As D is complex, Eqs|4^|48]and[32]require to be modified as follows 

^rr \ ^ ( Dr Di 

Wir Wii I '[ Di -Dr 


(44) 


(45) 


(46) 
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II 
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■uj 

+ IT^ 

ir 

■ Off ■ 

■u! 

+ IT^ 

-r rr 

■Go- 

Vk) 









-hIT^ ■ 

'* ir 

■ Gff 

■V,o 

(47) 

Fi -- 

II 

■ Gn 

■ uj 

+ IT- 

■ Gn 

■Uj 

+ K 

■Go- 

Vk) 










■ Gff 

■V,o 

(48) 


Notice that y,o e is not zero in this case. Both, the quadratic convex approximation and 
the analytical relaxed model can be extended for three-phase modeling. The problem increases 
in size but its solution is straightforward. 


5. Results 

Different simulations were performed in order to test the accuracy of the proposed method¬ 
ology. 
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The IEEE37 power distribution test system ll^ was used to evaluate the performance of the 
proposed methodology using a balanced equivalent modeling. Three generators were included 
in nodes 708, 732 and 744. Three optimization models were considered: Rectangular-OPE, 
quadratic programming formulation and relaxed analytical solution. Rectangular-OPE was im¬ 
plemented in GAMS 13111 as a non-linear/non-convex problem. The quadratic programming for¬ 
mulation was solved using the optimization toolbox of Matlab while the relaxed analytical formu¬ 
lation was solved by implementing Eq. dJSl l directly in Matlab. Results for S r{max) - S Umax) = 0.9 
are given in Table [T] Generated power for the quadratic programming formulation were exactly 
the same as for the relaxed analytical formulation (i.e inequality constrains were never active) 
on account of the high capacity of distributed generators. Moreover, this result was very close 
to the exact non-convex/non-linear programming formulation. Differences between the exact 
and the approximated model were less than 3%. Total active power loss of the system with¬ 
out distributed generation was Fi = 0.0298 while the optimized system was Fi = 0.0089. A 
considerable reduction. 


Table 1: Comparative results of the OPF on the IEEE37 test system with P^ax = Qmax = 0.9 


Variable 

Rectangular-OPE 

(GAMS) 

Quadratic convex 
model (MATLAB) 

Analytical relaxed 
model (MATLAB) 

§708 

0.849-i-j0.400 

0.8580-Hj0.4153 

0.8580-Hj0.4153 

§732 

0.042H-j0.028 

0.0420-Hj0.0210 

0.0420-Hj0.0210 

§744 

0.368H-j0.191 

0.388lH-j0.1913 

0.3881-Hj0.1913 

Pl 

0.00876 

0.0089 

0.0089 


A second simulation was performed on the same test system but this time the total capac¬ 
ity of generation was decreased to Fmax = Qmax - 0.7. Results are given in Table |2] Relaxed 
Analytical solution was not as good as in the previous case on account of the activation of in¬ 
equality constrains. Nevertheless, the quadratic programming formulation was still a very good 
approximation with an error of less than 3%. 


Table 2: Comparative results of the OPF on the IEEE37 test system with Pmax = Qmax = 0.7 


Variable 

Rectangular-OPE 

(GAMS) 

Quadratic convex 
model (MATLAB) 

§708 

0.700-Hj0.400 

0.700-Hj0.415 

§732 

0.147-bj0.028 

0.154-i-j0.021 

§744 

0.409-bj0.190 

0.413H-j0.191 

Pl 

0.0088 

0.00896 


The analytical relaxed model as initialization methodology for non-linear programming al¬ 
gorithms was also evaluated on the IEEE37 test system. Table [3] shows the number of iterations 
with and without initialization for different solvers available in GAMS. The number of iterations 
was drastically reduced in some cases (for example in IPOPT solver). 

Accuracy of the proposed approximation was analyzed in a more general context by a set of 
simulations over 1000 randomly generated radial distribution systems with parameters given in 
Table |4] Test system were compelled to be radial, not as a condition of the methodology, but as 
manner to simply the algorithm for random generation of test cases. 
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Table 3: Comparative results of the OFF on the IEEE37 test system using different solvers 


Solver 

With 

Initialization 

Without 

Initialization 

MINOS 

4 

6 

KNITRO 

6 

9 

CONOFT 

16 

17 

IFOFT 

9 

75 

LINDO 

18 

50 

FATHNLF 

1 

47 


Table 4: Parameters for the generation of random distribution systems 


Parameter 

Min 

Max 

Number of nodes (n) 

30 

60 

Line resistance 

0.001 

0.0170 

Line inductive impedance 

0.001 

0.0170 

Line capacitive admittance 

0.000 

0.0002 

Load power 

0.000 

0.2000 

Load power factor 

0.7 

1.0 

Constant power loads 

50%of« 


Distributed generators 

5%ofn 



First, Total power loss was calculated on each test system without distributed generation 
using a back-forward sweep algorithm 1^ . Then, an OFF was estimated using the proposed 
methodology. Next, a new power flow was calculated using the exact model. After all these 
calculations were performed on each randomly generated test system, a histograms were drawn 
to analyze the results. 

Figure |2] depicts the improvement in terms of power losses between the systems without 
distributed generation and the systems optimized by the proposed approximation of the OFF 
In most of the cases, improvement was more than 50%. As expected, distribution generation 
improves the losses. 

Accuracy of the solutions was demonstrated by an histogram of the error in power losses {ep) 
and voltages (ey) as shown in Figs[3]and|4] Fower error was less than 5% and voltage error was 
less than 2% in more than 80% of the cases. Notice that 2% of error does not mean 2% of voltage 
drop in the system. In these cases, the proposed approximation could be accepted as the solution 
of the OFF, whereas in the other cases it can be used as good initial approximation for a more 
accurate optimization model. 

High errors in other cases can be explained by their loading conditions. In those cases, the 
minimal voltage along the feeder is lower than 0.7 as shown in Fig |5] Fortunately, most of the 
power distribution systems (as well as the randomly generated test systems) have a good voltage 
profile and hence the proposed approximation is valid. 

Finally, the three-phase formulation of the optimal power flow was tested on the IEEE37 
unbalanced power distribution test system ifsol] . Balanced generation was imposed. Results 
are summarized in Table |5] Comparison in terms of maximum voltage error and losses error 
demonstrate the accuracy of the methodology even in the three-phase case. 
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Figure 2: Histogram the improvement in terms of power losses for 1000 randomly generated test power distribution 
systems 
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Figure 3: Histogram the error of power losses between the analytical relaxed methodology of OFF and the exact model 
for 1000 randomly generated test power distribution systems 


6. Conclusions 

A quadratic approximation of the OPF for power distribution systems was proposed. This 
approximation is based on the linear power flow for distribution systems. Convexity guarantees 
not only a unique solution but a fast, easy and efficient implementation in commercial optimiza¬ 
tion toolboxes. Results are very close to the Rectangular-OPF which is a non-linear/non-convex 
problem. 

An analytical solution for the unconstrained problem was also developed. This solution can 
be used as initial point for the non-linear programming formulation. Simulations using different 
commercial solvers available in GAMS shown the importance of a good initial point to accelerate 
convergence. 

Simulation results over a large set of random generated test distribution systems demonstrated 
the accuracy of the approximation. In most of the cases, the approximated solution had a very 
low error compared to the back-forward swing algorithm. 

The methodology was extended to three-phase systems. Additional constrains such balanced 
generation was strait forward included in the methodology. Y-connected loads as well as A- 
connected loads and single phase generations were included in the model. 

The use of a complex expansion instead of a Taylor series for linearization of the system, 
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Figure 4: Histogram the maximum voltage error between the analytical relaxed methodology of OFF and the exact model 
for 1000 randomly generated test power distribution systems 
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Figure 5: Histogram of the minimal voltage calculated by the back-forward sweep algorithm for 1000 randomly generated 
test power distribution systems 


although equivalent for most of practical purposes, can allow an easy an efficient modeling and 
implementation of the algorithms. 
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